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Abstract 

This paper is a continuation of [2] where a new model of biological invasions 
in the plane directed by a line was introduced. Here we include new features 
such as transport and reaction terms on the line. Their interaction with 
the pure diffusivity in the plane is quantified in terms of enhancement of 
the propagation speed. We establish conditions that determine whether the 
spreading speed exceeds the standard Fisher KPP invasion speed. These 
conditions involve the ratio of the diffusivities on the line and in the field, the 
transport term and the reactions. The asymptotic behavior for large diffusions 
or large transports is obtained. We also discuss the biological interpretation 
of these findings. 

1 Introduction 

This paper is a sequel of [2] in which we introduced a new model to describe biological 
invasions in the plane when a strong diffusion takes place on a line. The purpose 
of this model is to understand the effect of a line interacting with a homogeneous 
environment. This type of questions arises for instance in studying the propagation 
of diseases directed by roads [6] or the movements of animal populations in the 
presence of pathways allowing for a more rapid movement. An example of the latter 
is provided by recent observations of movements of wolves along seismic lines in 
Western Canada [5]. In [2] we derived the asymptotic speed of spreading in the 
direction of the line. There we showed that for low diffusion the line has no effect, 
whereas, past a threshold, the line enhances global diffusion in the plane and the 
propagation is directed by diffusion on the line. Moreover, the propagation velocity 
on the line increases indefinitely as the diffusivity on the line grows to infinity. 
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The goal of the present paper is to include new features in this model such as 
transport and reaction on the road and to understand the resulting new effects. 
These require new developments that turn out to make more transparent the case 
considered in [2] as well. 

Consider the line {(x, 0) : x G R} in the plane R 2 ; for the sake of simplicity, 
we will refer to the plane as "the field" and the line as "the road". For a single 
species, we consider a system that combines the density of this population in the 
field v(x,y,t) and the density on the line u(x,t). The main questions that we want 
to understand are the effects of a transport term qd x u as well as a decay rate p > 
on the road. An invasive species that can be carried by streams of water is an 
example of a situation where such additional terms are required. The transport and 
the decay rate are considered here to be uniform, that is, q and p are constant. The 
equations for u and v then read: 

d t u — Dd xx u + qd x u + pu = uv(x, 0, t) — pu x G R, t > 
d t v - dAv = f(v) (x,y)ett,t>0 (1.1) 

—dd y v(x, 0, t) = pu(x, t) — uv(x, 0, t) x G R, t > 0. 

Here, d, D, p, v > 0, p > 0, q E K and / G C 1 ([0, +oo)) satisfies the usual KPP type 
assumptions: 

/(0) = f(l) = 0, / > in (0, 1), / < in (1, +oo), f(s) < f'(0)s for s > 0. 

(1.2) 

Note that transport and pure decay only occur on the road and the question is to 
know how these interact with the diffusivity and growth in the field. We combine 
the system with the initial condition 

u\ t =o = u in M 

where u$, vo are always assumed to be nonnegative and continuous. 
Let ck denote the KPP spreading velocity [4] in the field: 

c K = 2y/df'(0). 

This is the asymptotic speed at which the population would spread in the open 
space - i.e., when no line is present (see [1]). 

We say that (1.1) admits the asymptotic speeds of spreading wf (in the direction 
ei = (1,0)) if for any solution (u,v) starting from a compactly supported initial 
datum (uo,vq) ^ (0,0), and for all e > 0, the following hold true: 



lim sup \(u(x, t), v (x, y, t))\ — 0, lim sup \(u(x, t), v(x, y, t))\ — 0, 

x<{w^ —e)t x>(w^. +e)t 

y>0 y>0 



Va>0, lim sup \(u(x, t), v{x, y, t)) - (U, V(y))\ = 0, 

( w t +e)t<x<(wt -e)t 
0<y<a 
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where (U,V) = (U,V(y)) is the unique positive, bounded, stationary solution of 
(1.1). 

Thus, the first step for proving the existence of the asymptotic speeds of spread- 
ing consists in deriving a Liouville-type result for stationary solutions. The existence 
of the asymptotic speeds of spreading implies in particular that 

Vc [w~, wt], lim (u(x + ct, t),v(x + ct, y, t)) = (0, 0), 

Vc E (w~,w+), lim (u(x + ct, t),v(x + ct, y, t)) = (U, V(y)), 

locally uniformly in (x, y) E Q. This kind of weaker formulation is sometimes used 
in the literature as the definition of the asymptotic speed. 

In [2], we proved that, in the case q,p — 0, (1.1) admits asymptotic speeds of 
spreading wf. We further identified a threshold situation: D/d = 2 below which 
wf = ±ck and above which = —w~ > ck- The first question we address in the 
present paper is how this threshold is modified by the presence of the transport and 
decay terms q and p. We derive the following 

Theorem 1.1. Under the assumption (1.2), problem (1.1) admits asymptotic speeds 



of spreading wt 1 (in the directions ±ei ). Moreover, if — < 2 + 

1 1 d ~ f (0) ^ v^f(0) 7 

^/ien = ±ck, else w+ > ck (resp. w~ < —ck)- 

This result lends itself to a biological interpretation which we now describe. 
The effect that we want to understand is to know when transport, diffusivity and 
mortality on the road combine to yield a spreading velocity larger than ck, the KPP 
invasion speed in the open field. For definiteness, let us consider propagation to the 
right (that is, in the direction ex). 

Let us first notice that the threshold for D/d after which the effect of the road 
is felt is raised by a decay rate p > and lowered by a transport q > on the road. 
On the contrary, when q < 0, the overall spreading to the right is hampered. In the 
absence of mortality on the road, that is when p — 0, the threshold condition reads 

D >2- * = 2 fl- q 



d ' ~ \fdW) V ck 
The above formula can also be written as 

D 



q>CK{1 '2d 

A remarkable consequence of this formula is that a transport q larger than ck is 
sufficient to enhance the overall propagation, no matter what the diffusivity ratio 
is. This explains for instance how a river may help the invasion of a parasite [3]. If 
q < ck(1 — D/2d) then the spreading is at the usual KPP velocity. This yields an 
interesting interpretation when considering the propagation against the direction of 
the transport, that is, when q < 0: for — q large enough, the asymptotic speed of 
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spreading in the direction +e 1 in the field is the KPP invasion speed. For instance, 
asymptotic propagation upstream against a river flow in the neighbouring field is 
unaffected by the river but downstream propagation, in the direction of the flow, 
can be significantly enhanced. 

We will actually carry out the whole study of system (1.1) in the case where the 
mortality term — pu is replaced by a more general reaction term g{u): 

d t u — Dd xx u + qd x u = uv(x, 0, t) — pu + g(u) t > 

d t v - dAv = f(v) (x, y)en, t>0 (1.3) 

-dd y v(x, 0, t) = fiu(x, t) — uv(x, 0, t) x e R, t > 0. 

The hypotheses are: /, g G C 1 ([0, +oo)) and satisfy 

/(0)=/(l) = 0, />0in(0,l), /<0in(l,+oo), (1.4) 

g(0) = 0, 3S> 0, g(S) < 0, (1.5) 

s i — y f(s)/s, s i — y g{s)/s are nonincreasing. (1-6) 

Condition (1.6) holds if / and g are concave. It implies that / and g are of KPP 
type: f(s) < f'(0)s and g(s) < g'(0)s. 

Theorem 1.2. Assume that (1.4)-(1.6) hold. Then: 

(i) (Liouville-type result). Problem (1.3) admits a unique, positive, bounded, sta- 
tionary solution (U, V). Moreover, U = constant and V = V(y). 

(ii) (Spreading). Problem (1.3) admits asymptotic speeds of spreading wf . 

(Hi) (Spreading velocity). If — < 2 — %— - =F — , % , then w^r = ±ck- Otherwise 

1 J v ; d ~ f(oy v^oj 

wt > c K (resp. w*_ < -c K ). 

Notice that, in Theorem 1.1, s \— > f(s)/s is not assumed to be nonincreasing. 
This hypothesis is only required in the proof of statement (i) of Theorem 1.2. An 
alternative hypothesis is given in Proposition 3.4 below. It is to be noted that the 
Liouville-type result may fail if s i— >■ f(s)/s is not nonincreasing. In this case, we 
still get a convergence result, but for initial data satisfying uq < vjp, vq < 1, and 
the convergence holds to the minimal positive, stationary solution (see Remark 3.5). 
Furthermore, the spreading speeds are still well defined and statements (ii) and (iii) 
of Theorem 1.2 hold true. 

When there is no transport on the road, there holds w~ = and the threshold 
condition given by Theorem 1.1 part (iii) becomes 

D < 2 g'(0) 

d - f(oy 

This allows us to understand the - somewhat mysterious - factor 2 in the threshold 
condition of [2]. Indeed, when g = 0, we recover the known condition D < 2d. We 
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further note that, if /'(0) = g'(0), the threshold condition becomes D = d - which 
is what one would have expected. Thus, when the same reaction occurs on the road 
and in the field, the diffusion on the road enhances the speed of propagation as soon 
as the diffusion on the road is larger than in the field. Lastly, from the above formula 
it follows that if g'(0) > 2/'(0) then the enhancement due to the road always occurs, 
whatever the diffusivity ratio is. 

In the last part of the paper, we investigate the limits of the asymptotic speeds 
of spreading wf as the diffusion and the transport on the road tends to +00. We 
derive the following 

Theorem 1.3. As functions of the variables D and q respectively, wf satisfy 



with h > independent of q and < k < 1 independent of D. 

The above result implies that, in the case of no reaction or of pure mortality on 
the road, a strong transport q allows the species to spread only at a portion k < 1 
of the speed of the flow. The expression of k is explicitly given in the proof of the 
theorem. 

The paper is organized as follows. In the next section, we start with recall- 
ing the results of [2] concerning the well-posedness of the Cauchy problems and 
the comparison principles. Then, we show that the large time limits of solutions 
emerging from non-trivial, bounded initial data are trapped between two positive, 
bounded, stationary solutions which do not depend on x. In Section 3, we derive the 
Liouville-type results for such class of solutions, first for (I.I) and then for (1.3). The 
large-time behaviour of solutions is thereby characterized, at least locally uniformly 
in space. In particular, Theorem 1.2 part (i) follows. 

In section 4, we investigate the spreading property of solutions, that is, the set 
where they converge to or to the positive, stationary solution. We give the proof 
of Theorem 1.2 parts (ii) and (iii) insisting on the main differences with [2]. 

Section 5 is dedicated to the proof of Theorem 1.3. 

2 Long time behaviour 

Throughout this section, we assume that (1.4), (1.5) hold. The unique solvability 
for the Cauchy problem associated with (1.3) in the class of nonnegative, bounded 
functions follows from the same arguments as in [2], Proposition 3.1. There, the case 
g = is treated, but a slight modification of the arguments allows one to handle the 
presence of the nonlinear term g. The existence result is obtained by first solving 
the Cauchy problem for u with v = and then the one for v with the so obtained 
u. Repeating this procedure, starting with the new v, one constructs an increasing, 
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bounded sequence (u n ,v n ) which eventually converges to a solution of (1.3). The 
uniqueness result is a consequence of the comparison principle between sub and 
supersolutions (i.e., pairs of functions satisfying the system with the "=" replaced 
by "<" and ">" respectively). We recall two versions of the comparison principle 
that will be used several time in the sequel: 

Proposition 2.1 ([2]). Let (u,v) and (u,v) be respectively a subsolution bounded 
from above and a supers olution bounded from below of (1.3) satisfying u <u and 
v <v at t = 0. Then, either u < u and v < v for all t, or there exists T > such 
that (u,v) = (u,v) for t < T. 

The second version deals with subsolutions in a generalised sense. 

Proposition 2.2 ([2]). Let E C R N x (0, +oo) and F C Q x (0, +oo) be two open 
sets and let (ui,vi), (u 2 ,v 2 ) be two subsolutions of (1.3) bounded from above and 



ui < u 2 on (dE) n (R N x (0, +00)), v x < v 2 on (dF) n (ft x (0, +00)). 
If the functions u, v defined by 



u(x, t) > u 2 (x, t) =>• v(x, 0, t) > V\(x, 0, t), 

v(x, 0, t) > v 2 (x, 0, t) u(x, t) > u±(x, t), 

then, any supersolution (u,v) of (1.3) bounded from below and such that u < u and 
v<v at t = 0, satisfies u<u and v<v for all t > 0. 

We now derive a result which gives a preliminary information about the long 
time behaviour of solutions. 

Lemma 2.3. Let (u,v) be the solution of (1.3) starting from a bounded initial datum 
(uo,vo) ^ (0,0). Then, there exist two positive, bounded, x -independent, stationary 
solutions (Ui,Vi) and (U 2 ,V 2 ) of (1.3) such that 



satisfying 





satisfy 



U\ < lim inf u(t, x) < limsup u(t, x) < U<< 



2- 



Vi(y) < lim inf v(t, x, y) < limsupw(t, x,y) < V 2 (y) 



locally uniformly in (x,y) G Q. 
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Proof. Let S be the constant in (1.5). The pair (U, V) defined by 



(U,V) 



max 



( 



^0 || oo + S ||Vo||oo + l 



V 



) 



is a supersolution of (1.3) which is larger than (u,v) at t — 0. Let (w, v) be the 
solution of (1.3) with initial datum (U, V). Using the comparison principle given by 
Proposition 2.1, we see that u and v are nonincreasing in t. Thus, (u, v) converges 
as t — > +oo to a stationary solution (U 2 , V 2 ) of (1.3) satisfying 



locally uniformly in (x, y) G f2. Furthermore, by translation invariance of the prob- 
lem in the x-direction, we see that (U 2 , V 2 ) does not depend on x. 

We now construct the pair (Ui, Vi). Take R > large enough in such a way that 
the principal eigenvalue of —A in B R C M 2 with Dirichlet boundary condition is less 
than /'(0)/(2d). The associated principal eigenfunction (p R satisfies 



Hence, for e > small enough, the function eipR satisfies — dA(eipR) < f(e(pn) in 
B R . We extend (f R to outside B R and we define V_(x,y) := E(p R (x,y — R — 1). 
Thus, (0,Y_) is a generalised subsolution of (1.3). The strong comparison principle 
given by 2.1 implies that u and v are positive for t > 0. Hence, up to decreasing e if 
need be, we have that (0, V_) is below (u,v) at, say, t — 1. Let (tt, w) be the solution 
of (1.3) starting from (0, V_) ait — 1. Using the comparison principle for generalised 
subsolutions - Proposition 2.2 - we see that u and v are nondecreasing in t. If one 
of them were not strictly increasing, the strong comparison principle of Proposition 
2.1 would imply that (u,v) is constant in time, which is impossible because (0,U) is 
not a solution of (1.3). Thus, as t — > +oo, (u,y) converges to a stationary solution 
(Ui,Vi) of (1.3) satisfying 

< U\{x) < liminf u(t, x), Y_{x,y) < V\(x,y) < liminf v(t, x, y), 

locally uniformly in (x, y) G f2. It remains to show that (U±, V\) does not depend on 
x. Since V_ is compactly supported, there exists k > such that (Ui, V\) is above the 
translated by any h G (—k, k) in the x-direction of (0, VI). By translation invariance 
of the problem, the solutions of (1.3) emerging from these initial data coincide with 
the translated by h G (—k, k) of (u,v). We then infer, by comparison, that (U\, Vi) 
is above the translated by h G (—k, k) in the x-direction of itself. Namely, it does 
not depend oni. □ 



limsup u(t, x) < U 2 (x), 



]imswpv(t,x,y) < V 2 (x,y), 



dA<p R <-f(0)<p R mB R . 
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3 Liouville-type result for 1-dimensional solutions 



In this section, we derive a Liouville-type result for stationary solutions of (1.3) 
which do not depend on x. Namely, we will show that the problem 

'U = constant, V = V(y) 
< -dV" = f(V), y>0 
vV(0)=fiU-g(U) 
-dV'(0) = g(U). 

admits a unique positive, bounded solution. This will imply the general Liouville- 
type result, Theorem 1.2 part (i), because, by Lemma 2.3, any positive, bounded, 
stationary solution of (1.3) lies between two positive, bounded solutions of (3.1). 

We start with considering the pure mortality case g(U) = —pU. The proof is 
much simpler in this case. Problem (3.1) reduces to 

-dV" = f(V), y > (3.2) 
^'(0) = ^-V(0). 

Proposition 3.1. Under the assumption (1.4), problem (3.2) admits a unique pos- 
itive, bounded solution. 

Proof. Let V be a positive, bounded solution of (3.2). It is straightforward to check 
that V necessarily satisfies < V < 1, V(+oo) = 1 and V'(+oo) = 0. Thus, 
multiplying the second equation of (3.2) by V and integrating by parts between 
and +oo we get 

Examining the function 9 defined by 
we see that 

0(0) <0, W > 0, 0'>Oin(O,l). 

There exists then a unique value do G (0, 1] such that 0(ao) = 0. Hence, V(0) = uq. 
This proves the uniqueness of positive, bounded solutions of (3.2). It is also easy to 
verify that the solution V of (3.2) with initial datum V(0) = cr is actually positive 
and bounded. Assume indeed by contradiction that this is not the case. One can 
then find f > such that either V(£) = 1 and V'(£) > 0, or V(£) < 1 and V'(£) = 0. 
Owing to (1.4), both cases are ruled out by the following equality: 

[ Vi ° f(s)ds = d -\{V\U)f - (V"(0) 2 ] = f f(s)ds - d -{V\i)f. 

□ 
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This concludes the proof of the Liouville-type result for stationary solutions of 
(1.1). 

We now pass to the case of a general reaction term on the road. 

Theorem 3.2. Under the assumptions (1.4)-(1.6), problem (3.1) admits a unique 
positive, bounded solution. 

The existence result is contained in Lemma 2.3, that we proved using a sub and 
supersolution argument. Let us present a more explicit construction inspired by the 
shooting method. 

Proof of the existence part of Theorem 3.2. For U EM., let Vy denote the associated 
solution of (3.1). Consider the following set: 

U := {U > : Vy > 0, Vy(y) > 0}. 

This set is nonempty because, for U > max(z///i, S) (S being the constant in (1.5)) 
the function Vy satisfies Vu(0) > 1, V^(0) > 0. It follows that Vy is nondecreasing 
and then positive. We then define 

U* := MU. 

Suppose by contradiction that U* = 0. We can then take U EU close enough to 
in such a way that Vu(0) < 1. Notice that V v (0) > 0, because otherwise Vy would 
not be positive. Call r\ the first point where Vy reaches the value 1, if it exists, 
else set i] := +oo. In the second case, the function Vy lies in (0, 1) and then, since 
dV{l = —f(Vy) < 0, it satisfies Vy(+oo) = 1, V^(+oo) = 0. Hence, in both cases, 
V{j{r}) > 0. We derive 



1 P 1 P V hi) 

vm>- d f /«,<»))*> 5 jf ^mwy 



dV{,(0) 



Choosing U EU small enough then leads to a contradiction, because Vy(0), V^(0) — > 
as U ->• 0. Therefore, U* > 0. 

We claim that Vy* is positive and bounded. The continuity of solutions of the 
Cauchy problem with respect to the initial datum yields Vy* > 0. Assume by way 
of contradiction that there exists y$ > such that Vy*(yo) = 0. We necessarily 
have that V u *(y ) = 0, because otherwise Vy* would be negative somewhere. Hence, 
Vy* = by the uniqueness of solutions of the Cauchy problem. But then U* = 0, 
which is impossible. If Vy*(0) = then 

V£.(0) = -g(U*)/d = -fiU*/d < 0, 

which is again a contradiction. We have shown that Vy* > on [0, +oo). 

Assume now by contradiction that Vy* is not bounded from above. There exists 
then y > such that 

Vy*(y ) > 1, V u ,(y o )>0. 
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Since V v — > V v * as U — >■ £/* in C/ oc ([0, +00)), for [/ close enough to U* we have that 
min V v >\ min > 0, Vu(y ) > 1, V{j{y ) > 0. 

[0,2/ ] ^ [0,2/o] 

It follows that Vjy is increasing in [y , +00) and then min Vjj > 0. This means that 
U 6W, which contradicts the definition of U* . Therefore, (U*,Vu*) is a positive, 
bounded solution of (3.1). □ 

Proof of the uniqueness part of Theorem 3.2. Assume by way of contradiction that 
(3.1) admits two distinct positive, bounded solutions (U 1: Vi) and (U 2 ,V 2 ). The 
uniqueness of solutions of the Cauchy problem yields U\ 7^ U 2 ; say, U\ < U 2 . Since 
(1.6) implies that the function G : R+ — > R defined by G(s) :— ji — g{s)/s is 
nondecreasing, we obtain 

i/Vi(0) = G(C/i)C/i < G(C/ 2 )C/i < ^2(0), 

and equality holds if and only if V 2 (0) = 0. But if V 2 {0) = then -dV 2 \0) = fiU 2 > 
and then V 2 would not be nonnegative. Thus, Vi(0) < V^O). We argue differently 
depending on the fact that V\ = V 2 somewhere or not. 

Case 1) Vi(y) = V 2 {y) for some y > 0. 
Let 7] be the smallest y > such that V\{y) = V 2 (y). Namely, V\ < V 2 in (0,r?) and 
V\{rj) = V 2 (r]). By the uniqueness of solutions of the Cauchy problem we infer that 
V{(r]) > V^irj). Multiplying by V 2 the equation satisfied by V" and by V\ the one 
satisfied by V 2 " and integrating over (0,77), we get 

<V((0)V 2 (0)- VZ(0)Vi(0). 
The first integral above is nonnegative by (1.6). Thus, 

Vi'(0)y 2 (0) > ^'(0)^(0). (3.3) 

This inequality reads 

giUJifJiUt - g(U 2 )) < giUJifiU, - g^)), 

that is, 

g{Ui)U 2 < g^Ut. 

This contradicts (1.6) because U\ < U 2 . 

Case 2) V^y) < V 2 (y) for all y > 0. 
Notice that V\ and V 2 cannot attain the value 1 with a nonzero first derivative, 
because they would be either unbounded or negative somewhere. It follows from 
the uniqueness of solutions of the Cauchy problem that they are either identically 
equal to 1, or below 1, concave and increasing or above 1, and decreasing. In any 
case they satisfy 

Vi G {1, 2}, hm V t (y) = 1, lim V?(y) = 0. 
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Thus, the same computation as in the case 1 yields 

= Vl(0)V 2 (0)-Vi(0)V 1 (0). 

We know from (1.6) that the function s i— >■ f(s)/s is nonincreasing. Moreover, since 
it is positive in (0, 1) and negative in (1, +oo), it cannot be constant in a left or in 
a right neighborhood of 1. It follows that f(Vi)/Vi > /(V^/V^ on [0,+oo), with 
strict inequality somewhere. As a consequence, the first integral above is strictly 
positive. That is, (3.3) holds. We then get a contradiction by arguing as in the 
previous case. □ 

We now derive some properties of positive, bounded solutions of (3.1). 

Proposition 3.3. Assume that (1.4), (1.5) hold and that s h-> g(s)/s is nonincreas- 
ing. Let 5* := inf-f^ > : g(S) < 0}. Then, any positive, bounded solution (U, V) 
of (3.1) satisfies 

~ A* [V<1, ~ A 4 \^>1. 

Proof. Since — cH/" = f(V), V cannot attain the value 1 without being constant. 
The following implications are then easily obtained: 

V = 1 on [0, +oo) V(0) = 1 <£> V'(0) = g{U) = S* < U =-, 

At 

V < 1 on [0, +oo) V(0) < 1 V'(0) > ^ #([/) < S* < U <-, 

fi 

V > 1 on [0, +oo) V(0) > 1 ^ V'(0) < ^ #([/) > - < U < S*. 

At 

The result follows. □ 

It is not hard to construct examples with s h- > f(s)/s not nonincreasing for which 
the Liouville-type result fails. However, as shown by Proposition 3.1, this hypothesis 
can be dropped if g satisfies some suitable conditions. The weaker sufficient condition 
we are able to obtain is expressed in terms of the following quantity: 

Sm '■= min{s > : g' < in [s, +oo)}. 

Since Sm = when g(s) = —ps, the following result generalizes Proposition 3.1. 

Proposition 3.4. Assume that (1.4), (1.5) hold, that s >->■ g(s)/s is nonincreasing 
and that 

u + g(S M ) 

S ■ 

At 

Then, problem (3.1) admits a unique positive, bounded solution. 
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Proof. Let (U, V) be a positive, bounded solution of (3.1). We know that V^+oo) = 
1 and V'(+oo) = 0. The same integration by parts as in the proof of Proposition 
3.1 yields 9(U) = 0, with 

9(a) := g 2 (a) + 2d f(s)ds, G(a) := (^a - g(a))/u. 

Thus, U satisfies 



((G(U) 
\9(U) = 



' ^ * ° (3.4) 
0. V 1 



It is not hard to show that, conversely, if U > satisfies (3.4) then the associ- 
ated solution of (3.1) is positive and bounded. The function G satisfies G(0) = 0, 
G(+oo) = +oo and it is strictly increasing in the set where it is positive. Let 
Si > be such that G(Si) = 1. Set S* := inf{S > : g(S) < 0} and then call 
a := min(Si, S*), b := max(Si, S*). For U > b, G{U) > 1 and g(U) < 0. Hence 
U > b satisfies the first condition in (3.4) only if g(U) = 0, but then 6{U) < 
because / < in (1, +oo). If < U < a then (G(U) - l)g(U) < 0. This shows that 
(3.4) has no solution outside [a, b}. Since 6(S*) < and 9 (Si) > 0, the function 9 
vanishes somewhere on [a, b]. Moreover, (G — l)g > there. Therefore, (3.4) admits 
solution on [a, b]. We conclude the proof by showing that 9 is strictly monotone on 
[a,b]. If S* < Si then g 2 is nondecreasing and G is strictly increasing and smaller 
than 1 on [a, b], whence 9 is strictly increasing. Consider now the case S* > Si. The 
function G is strictly increasing and larger than 1 on [S^S*]. On the other hand, 
by hypothesis, 

G(S M ) = ^ Sm ~J {Sm) < 1. 

Whence, Sm < and thus g 2 is decreasing on [Si, S*]. We eventually infer that 9 
is strictly decreasing on [Si, S*]. □ 

Remark 3.5. Under the assumptions of Proposition 3.3, if (3.1) admits multiple 
positive, bounded solutions, then there exists a minimal one (U_,V_) among them. 
Moreover, (U_,Y_) satisfies U_ > uf/i, V_ > 1 and it attracts solutions of (1.3) starting 
below (v/n, 1). Indeed, if non-uniqueness occurs, then Proposition 3.4 yields 

S, > s M > > 

Whence any positive, bounded solution (U, V) of (3.1) satisfies U > u//j, and V > 1 
due to Proposition 3.3. Therefore, by comparison, if (u,v) is a solution of (1.3) with 
an initial datum (uo,vo) below (v/n, 1), then 

limsup u(t, x) < U, limsup v(t, x, y) < V(y), 

locally uniformly in (x,y) G Q. On the other hand, if (uo,vo) ^ (0,0), Lemma 2.3 
implies that equality holds in the above expressions for one of these solutions (U, V). 
This is the minimal solution. 
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4 Propagation 



In this section we prove statements (ii) and (iii) of Theorem 1.2. We concentrate 
on propagation to the right, since propagation to the left is obtained by replacing 
q with —q. The general plan of the proof is that of [2]. First, we look for plane 
waves (or exponential solutions) for the linearised system around (0, 0). These have 
a critical velocity w*. It then follows from the KPP assumption that solutions of 
the nonlinear system spread at most with velocity w*. To prove that they spread 
at least with velocity w*, we look at the problem in a wide strip R x (0, L) and we 
construct compactly supported sub-solutions using again exponential solutions, but 
this time with complex exponents. Finally, we establish the condition under which 
the critical wave velocity is higher than the KPP speed. We do not repeat here all 
the arguments of [2], but we give details only of the points that are different. 

4.1 Exponential solutions 

Exponential solutions of the linearised system are looked under the form 

(u(t,x),v(t,x,y))^e- a ^- ct \l, 1 e-^), (4.1) 

with c > 0, 7 > and a, (3 G R (not necessarily positive). Namely, we look for c, a, 
[3 and 7 satisfying 

ca — Da 2 — qa — 1/7 — fj, + g'(0) 
ca - d(a 2 + (3 2 ) = /'(0) 
d^y/3 — fi — vy. 

The third equation yields 7 = /i/(z/ + d/3), which is positive iff /3 > —vjd. This will 
be assumed without further reference. Substituting we get 

ca-d(a 2 + (3 2 ) = f(O). 
The first equation of (4.2) in the unknown a has the roots 

P) = ^ (c - q ± V(c-g) 2 + 4D( X (/3W(0))) • 
where the function x is defined on (—v/d, +00) by 

v + dp 

Since \ is strictly increasing and tends to /i at +00, we see that ar^(c, (3) are real iff 

(c -q) 2 > AD(g'(0) -//), (4.3) 
(3 > (3(c) := x -\g>(0)-(c-q) 2 /4D). 
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Therefore, if c satisfies (4.3), the set of real solutions ((3, a) of the first equation of 
(4.2) is given by £(c) = £~(c) U £+(c), with 



E^c) :={(/?, a±(c,/3)) : > £(c)}. 

This is a smooth curve with leftmost point (/3(c), (c — q)/2D). The second equation 
in (4.2) admits real solutions iff c > ck, where ck '■= 2a/g(/ 7 (0) is the classical KPP 
critical speed for the second equation in (1.3). In the ((3, a) plane, it represents the 
circle T(c) of centre (0, c/2d) and radius 

r(c) = "^d— 

Let «S(c) denote the closed set bounded from below by £~(c) and from above by 
S + (c) and let Q(c) denote the closed disc with boundary T(c). Exponential functions 
of the type (4.1) are supersolutions of the linearisation of (1.3) iff (/3,a) G S(c)r\Q(c). 
It is easy to check that, for ck < c < c', £/(c) C £?(c'). Moreover, as c — > +oo, the 
radius r(c) of the disc C?(c) diverges and its bottom (0, (c— \/ c 2 — c 2 K )/2d) tends to 
(0,0). Consequently 

[J 0(c) = Kx (0,+oo). 

C>C K 

On the other hand, the sets S(c) are not increasing with respect to c, but the sets 
S(c) fl(Ix R + ) are. Indeed, under condition (4.3), a^(c, /3) > iff c > g or c < g 
and x(/3) > ^'(O)) an d m both cases 

2Dd c at(c, 13) = 1 + °~ q = > 0, 

y/ic-qF + WixW-m) ~ 

whereas a~^(c, (3) > iff c > q and < #'(0), which yields 

2Dd c a-Jc, p) = \ C ~ q < 0. 

y/ic-qy+wixm-m) ~ 

It follows that there exists > ck such that S(c) C\Q(c) ^ iff c > w*. Moreover, 
= ck iff (0, CK/2d) G S(c K ). Otherwise S{w*) D reduces to a point in 

(0, +oo) 2 , denoted by (/?*, a*). These notations will be kept in the following sections. 



4.2 Spreading 

Proof of Theorem 1.2 part (ii). Let (w,f) be the solution of (1.3) starting from a 
nonnegative, compactly supported initial datum (u ,v ) ^ (0,0). We prove sepa- 
rately that (u, v) spreads (towards right) at most and then at least at the critical 
velocity of plane waves w*. 

Step 1. (u,v) spreads at most with velocity w*. 
The definition of implies the existence of an exponential supersolution of the 
linearization of (1.3) (thus a supersolution of (1.3) by (1.6)) of the type (4.1), with 
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Figure 2: Case f£ G S(cjr). 



c = if*, 7, a > and /3 > 0. Since, up to translation in the direction e±, this 
supersolution is above (u,v) at time t = 0, the comparison principle yields 

lim sup \(u(x, t), v(x, y, t))\ = 0. 

t^ + OO a!> (,„ t+e )t 
V>0 

Step 2. (w, u) spreads at least with velocity 
We need to show that 

V£:,a>0, lim sup \(u(x,t),v(x,y,t) — (U,V(y))\ — 0, 

t^+CC 0<x<(w*-s)t 
0<y<a 

where (U,V) is the unique positive, bounded, stationary solution of (1.3). We first 
show that 

We > 0, 3c G - e, lim (u(x + ct, t),v(x + ct, y, t)) = (U, V(y)), (4.4) 

t— >+oo 

locally uniformly in (x, y) G Q. Then we will conclude applying the following lemma 
with c~ = 0, c + = c. 

Lemma 4.1. Let c~ < c + be such that any nonnegative, bounded, solution (u,v) ^ 
(0,0) of (1.3) satisfies 

lim (u(x + c ± t,t),v(x + t*t,y,t)) = (U,V(y)), 

t— >+oo 
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locally uniformly in (x,y) E fl, with U, V > 0. Then, such solutions satisfy 
Va > 0, lim sup \(u(x, t), v(x, y, t) - (U, V{y))\ = 0. 

*^ + °° c -t<*<c+t 
0<y<a 

In order to prove (4.4) it is sufficient to derive the analogue of Lemma 2.3 for the 
problem in the frame moving with speed c in the direction ei, for c arbitrarily close 
to w*. Then one concludes using Theorem 3.2 (or Proposition 3.1 in the case of (1.1) 
under the only assumption (1.2)). The only argument in the proof of Lemma 2.3 
which is modified by the change of frame is the existence of a compactly supported, 
stationary, strict subsolution (U_,V) for the linearised problem. If w* = ck then 
the construction of a function V_ such that (0, V_) satisfies the desired properties is 
standard (see Lemma 6.2 of [2]). In the case > ck, the bifurcation analysis of 
Lemma 6.1 of [2] applies and provides the subsolution (U_,V_) for c close enough to 
u>*. Since this is the core of the argument, for the sake of completeness we sketch 
below the proof of that lemma. 

One starts with penalizing f'(0) and g'(0) and considering the problem in the 
strip R x (0, L), with a Dirichlet condition on R x {L}: 



—Dd xx u + (q — c)d x u = vv(x, 0, t) + (g'(0) — 5 + /jl)u x E R 

-dAv - cd x v = (/' (0) - 5)v (x, y)ERx (0, L) 

v(x,L)=0 xER 

—dd y v(x, 0, t) = jiu{x) — vv(x, 0,t) x E R. 



(4.5) 



The presence of 5 « 1 can be viewed as a perturbation of the terms f'{0), g'(0); 
thus we can continue the discussion with 5 = 0. Exponential solutions are sought 
for in the form 

e- ax (l,7e-^ + 7e^). 

They exist for (/3,a) E T(c) fl Sl(c), where T(c) is the same as in the previous 
section and £z,(c) ^ s ^ ne union of two curves S^(c) parametrized by a = a^(c,/3). 
The functions have the same monotonicity in c as the functions cc 11 defining 
£(c); moreover, as L — > oo, — > locally uniformly in (3 > 0, together with 
their derivatives. Calling S L (c) the set between S£(c) and Sj(c), it follows that 
S L (c) fl Q(c) 7^ iff c > with — y as L — > oo. In particular, > c K for 
L large enough (and 5 small enough). By using Rouche's theorem one eventually 
finds, for c < close enough to w^, two values a, (3 E C\R such that the associated 
exponential function satisfies (4.5). Using its real part, one eventually obtain the 
compactly supported subsolution (U_,V_) for c arbitrarily close to w*. □ 

Proof of Lemma J^.l. Let (u,v) be as in the statement of the lemma and fix a > 
0, < e < U. Consider the solution (u l ,v l ) emerging from the initial datum 
(sup u, sup v), and a solution (u 2 ,v 2 ) ^ (0,0) with an initial datum (u^v 2 ) satisfying 

0< Mo <£/-£, SUpptiQ C [-1, 1], Wq = 0. 
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By hypothesis, for a > 0, there exists T > such that 

Vt>T, 0<y<a, |(^(c~t, t), v\c~t, y, t)) - (U, V(y))\ < e. (4.6) 

Set k := max(l, |c _ |T) and let T' > be such that 

Vt>T', \x\ < k, < y < a, \(u(x + c + t,t),v{x + c + t,y,t)) - (U,V(y))\ < e. 

For 1/2 < A < 1, call c := (1 - A)c~ + Ac + . Let s > 2T'. If (1 - A)s < T, writing 
cs = x + c + t with | a; | < k and t > T", we find that 

Vy e [0,a], |Kct,0,^,y,t))-(C/,y(y))| <£. 

Consider the case (1 — A)s > T. For (:r,y) G fi, we see that 

ul(x) < (U — < -u(a; + c + As, As) < -Up, v o ^ ^(^ + c + As, y, As) < t>o, 

and then the comparison principle yields 

u 2 (x + (1 - A)fTs, (1 - A)s) < u(x + cs, s) < u^a; + (1 - A)c~s, (1 - A)s) 

v 2 (x + (1 — A)c~s, y, (1 — A)s) < v(x + cs, y, s) < v (x + (1 — A)c~s, y, (1 — A)s). 

Whence, by (4.6), 

Vye [0,a], \(u(c-s,s),v(cS,y,s)) - (U,V(y))\ <e. 

We have therefore shown that (u, v) converges to (U, V) as t — > H-oo uniformly in 
the set (c~ + c + )t/2 < x < c + t, < y < a. The proof is completed by exchanging 
the roles of c~ and c + . □ 



4.3 Spreading velocity 

Proof of Theorem 1.2 part (in). We have shown in the previous section that the 
asymptotic speed of spreading coincides with the critical velocity of plane waves u>*. 
With the same notation as in Section 4.1, we know that w* = ck iff (0,ck/2cI) E 
S(ck), that is, iff a^(cK, 0) are real and satisfy 



a D (c K ,0) <7^< a^(c K ,0). 



(4.7) 



Requiring (4.7) is equivalent to 
ck c k - q 



2d 



which, in turn, rewrites 



2D 



<^V(c K -q) 2 -4Dg>(0), 



( Cx (7 -1 ) +g ) -( CK - q)2 - AD9 '( 0) - 

Notice that this condition automatically implies that a^(cK,0) are real. Whence, 
recalling that ck = 2^/df'(0), we find that (4.7) holds iff 



4^/'(0) + 4^/^ q - 8/'(0) < -V(0). 



This concludes the proof. 



□ 
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5 The large diffusion and transport limits of wf 



We now prove Theorem 1.3. As before, we focus on the speed towards right, recalling 
that it coincides with the critical wave speed w* defined in Section 4.1. 

5.1 Large diffusion 

Replacing c with c = y/Dc and a with a/y/D reduces (4.2) to the following system: 

-or + c = a = — + g (0) 

V VBJ v + dP (5.1) 

c Dc 

The first equation is satisfied for (/3,a) £ £(c — q/y/D), where £ = £~ U £ + is 
the curve defined in Section 4.1 with D = 1 and q — 0. The second equation is 
that of an ellipse Ed(c) in the (/3, a) plane, which is above the parabola P(c) of 
equation a = (/'(0) +<i/3 2 )/c. For c close to 0, S(c) is below P(c). Then, increasing 
c, S + (c) moves upward while P(c) moves downward and tends to the j3 axis as 
c — > +oo. There exists then a positive value h such that E(c) H P{c) ^ if and 
only if c > /i. Since Ti{c — q/yfD) and Ed(c) converge locally uniformly to E(c) and 
P(c) respectively as £> — >■ oo, and E is bounded in the a direction, it follows that if 
c < h then (5.1) has no solution for D large enough, whereas if c > h then (5.1) has 
solution for D large enough. Reverting to the original variable, we eventually infer 
that 




(a) c < h (b) c = h (c) c> h 

Figure 3: The asymptotics of as D — > oo. 

5.2 Large transport 

We consider now as a function of q. We know that = for — q large enough. 
Let us investigate the behaviour as q — > +oo. Notice first that, as q — > +oo, 
aj(c, /3) — >■ locally uniformly in c and uniformly in j3 > 0, from which we deduce 
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that — > +00. Let us take c = nq with k > and q — >■ +00. The set 
intersect if and only if there exists (5 > such that 



-L(, (K - 1) + ^ (K _ 1)2 + 4JW) _ g , (0)) ) > gzjgffig . 

If ac > 1, this inequality holds for any given (3 > 0, provided q is large enough. This 
shows that lim sup _>. +00 w*/q < 1. If/c < 1, the inequality implies that, as g — >• +00, 

x(/3)-</(0) m . 2f(0) + 2d/3 2 

+ o(l) > 



1-k k + v//t 2 - (4 + 4rf 2 /3 2 )g- 2 ' 

In particular, /3 cannot diverge as q — > +00, whence 

xW) - m + o(1) > m+jg 

1 — K _ K 

It follows that, if g'(0) > /i and k < 1, r(/tg) flS(/tg) = for g large enough. Instead, 
if g'(0) < /x, there is a threshold value 

such that r(/tg) fl S(/tg) is empty if k < and g is large enough and is nonempty if 
k> k and q is large enough. Theorem 1.3 follows. 

Remark 5.1. The computations in both the large diffusion and transport cases 
show that the scale at which the limit of w* as D, q — > +00 is affected by both terms 
is D ~ q 2 . This was of course expected by dimensional considerations. 
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